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The cold dark matter may be in a meta-stable state and decays to other particles with a very long 
lifetime. If the decaying products of the dark matter are weakly interacting, e.g. neutrinos, then it 
would have little impact on astrophysical processes and is therefore difficult to observe. However, 
such a decay would affect the expansion history of the Universe because of the change of the equation 
of state. We utilize a high-quality type la supernovae (SN la) data set selected from several resent 
observations and the position of the first peak of the Cosmic Microwave Background (CMB) angular 
spectrum given by the WMAP three-year data to constrain the dark matter decay-to-neutrino rate 
r = aF^, where a is the fraction of the rest mass which gets converted to neutrinos, and is the 
decay width. We find that > 0.7 x 10^ Gyr at 95.5% confidence level. 



I. INTRODUCTION 

The nature of dark matter is a great mystery of the 
Universe. Although it constitutes about a quarter of 
the total cosmic density, its many properties are still 
not known. Various kinds of dark matter candidates 
have been proposed, for instance, the weakly interact- 
ing massive particle (WIMP) is one of the most pop- 
ular candidates l]. However, the standard Cold Dark 
Matter (CDM) still have some problems at small scales, 
namely the halo cuspy profile problem and the missing 
satellite problem 0, IH, which motivated researchers to 
study other possibilities 

Decaying cold dark matter model (DCDM) is one of 
such possibility. It has been suggested that DCDM may 
help resolve the small scale problems of the standard 
CDM Furthermore, from the particle physics per- 
spective, it is also natural to think that the dark matter 
particle may be a meta-stable particle, which could de- 
cay into more stable particles. Such possibility could be 
realized if the dark matter is not the true Lightest Su- 
persymmetric Particle (LSP) [ol, or if the R-parity of the 
LSP is not exactly conserved [7|- 

If the decaying products are photons, electrons, quarks 
or other particles which are involved in the strong or 
electromagnetic interactions, then it is relatively easy to 
be observed. Stringent limit on radiative particle decay 
have been derived, e.g., from the reionization history of 
the Universe However, if the decaying products 

include only weakly interacting particles, e.g. neutrinos 
or other WIMPs, then it is more difficult to obtain such 
constraint. 

We consider a generic unstable cold dark matter par- 
ticle which decays invisibly, producing only weakly in- 
teracting particles. Within the standard model of parti- 
cle physics, the only such weakly interacting particle is 
neutrino. The unstable decaying particle may also pro- 
duce another weakly interacting particle which is beyond 
the standard model. For example, if the sneutrino is a 
meta-stable particle, it may decay to a neutralino and a 
neutrino, conserving both lepton number and R-parity. 

Since neutrinos have very small masses, those produced 



in such decays will be highly relativistic. We assume here 
that the recoil on the heavier decaying product is small 
so that it could still be treated as part of CDM. As the 
non-relativistic DCDM decays to relativistic neutrinos, 
the equation of state changes from to 1/3, and this will 
affect the expansion rate H . Observational constraint on 
relativistic component were discussed previously, see e.g. 
Ref. Tl, 12, 13]. 

Here, we constrain the invisible decay of dark mat- 
ter using two observations which probe the cosmic ex- 
pansion rate. The observations we use include (1) type 
la supernovae (SN la) data selected from several re- 
sent observations (involving Gold06 [l3|, SNLS [lB| and 
ESSENCE [HI); (2) the position of the first peak of 
the Cosmic Microwave Background (CMB) angular spec- 
trum £i given by the three-year WMAP observations 
[TtI [T§ |. These two observations cover a large redshift 
range (^SNaoi ~* ^■'^ ^^'^ -^cmb —>- 1048) and could give a 
tight and reliable constraint on the DCDM. Throughout 
the paper, we assume the dark energy is a cosmological 
constant and the geometry of the Universe is flat. 



II. CONSTRAINTS 

The evolution equations for the energy density of 
DCDM particle and its decay products p^, are given 

by 



Px + 

Pr + 411 Pr 
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where F ~ oT^, F^ is the decay width of the particle, 
< a < 1 is the fraction of the rest mass of the dark 
matter particle which goes to relativistic neutrinos during 
each decay. The Hubble expansion rate H is given by 



H' = 
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here pb is the baryon energy density which evolves as 
(1 -I- z)^, Pr is the relativistic particle energy density, 
which in the absence of decaying contribution would 
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evolves as (1 + z)**, and pA is the energy density of cos- 
mological constant. Substitute ^ = —H{1 + z)j^ into 
Eq. ([3]) , the redshift evolution of the energy density for 
different components can be solved numerically. We de- 
note flio — Pio/Pco, where pco is the critical density at 
z — 0, and take the radiation background (photon and 
neutrinos produced in the early Universe) in the absence 
of decaying contribution as Clro — 4.183 x lO'^h^^ where 
ho = i/o/lOO kms^ Mpc^ is the reduced Hubble con- 
stant. The baryon density is given by the CMB [13, UB] 
and big bang nucleosynthesis measurements [l9l | , though 
the actual value does not significantly affect our result. 
The cosmological model is then described by the param- 
eter set 

^ — (^xOi ^1^0, r, ho). 

The type la supernovae can be used as "standard can- 
dles". We have selected a sample of 182 high quahty 
SN la data from several resent observations, including 30 
HST supernovae, 47 SNLS supernovae from the Gold06 
data set, 60 ESSENCE supernovae, and 45 nearby super- 
novae from WV07 20]. The MLCS2k2 hght-curve fitter 
[2ll . [23 | is used to process this data set. The for SN 
la data is 
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where p,obs and ai are the observational modulus and its 
error, and fith is the theoretical modulus given by 

pif,j(z) = 51ogiodL(2;) + 25, (5) 

where the luminosity distance can be computed as 



dL{z\e) = (i + z) 



cdz' 
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In order to constrain the evolution at high redshift, we 
also use the position of the first peak of the CMB angular 
power spectrum li as measured by the WMAP three-year 
data. An alternative way of using CMB to constrain 
cosmic expansion history without invoking a full CMB 
calculation is to use the so called shift parameter R [2^ , 
which is a geometric test which assumes that the size 
of the sound horizon at recombination is fixed [24, 25]. 
However, the shift parameter was derived for assuming 
constant CDM density, which might not apply in our case 

The angular scale of the sound horizon at last scatter- 
ing, ^A = T^r{zis)/rs{zis), is the ratio of the comoving 
angular diameter distance of the last scattering surface 
r{zis) to the comoving size of the acoustic horizon rs{zis) 
at decoupling, which denotes both the size of sound hori- 
zon and the geometrical property of the Universe. Relat- 
ing to the cosmic expansion rate, this acoustic scale can 
be written as 5a 



^is cdz' 
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l/(l+^i.) 



(7) 



where zu is the redshift of last scattering, a is the scale 
factor and 
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-1/2 
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is the sound speed 



Here, fi-yo — il.ro/1-681 is the 



energy density of photon today. The value of Ia is re- 
lated to ii [23] , the first peak in the CMB angular power 
spectrum. Following Hu et al. [28l ] , we use the relation 



^1 =£^(l-0i). 



(9) 



where (pi is a phase factor, and if we set the spectral 
index n = 1 an 
fitting formula: 



index n = 1 and flboho = 0.02, then it is given by the 



6,= 0.267(^)°\ 



(10) 



Here = p^{zis) / pm{zis) and Pm{zis) = Pb{zis) + 
Pxizis)- This fitting formula is consistent with the other 
forms [29| and is a good approximation for our model, 
especially for the long lifetime DCDM. 

According to the measurement of the three-year 
WMAP data [ll[l3, we have £i = 220.7±0.7 and finally 



220.7)2 



0.72 



(11) 



For the combined SN la and CMB data set, the is 
given by 



X 



(12) 

We employ the Markov Chain Monte Carlo (MCMC) 
technique to simulate the probability distributions of the 
parameters in the model. The following priors on the 
parameters are adopted: il^o & (0,1), ^1^0 & (0,1), 



(6) logio(r/10- 



e (-7.5,6.5) and ho £ (0.2,0.9). We 



also assume that the energy density of each components 
is positive, fl^o £ (0, 1 — f^^o ^ ^bo — ^ro) is set so that 
Ha = 1 — fly^o — ^uo — ^bo — ^ro > 0. We generate 
six MCMC chains, and each chain contains about one 
hundred thousand points after convergence is reached. 
Finally, these chains are thinned to produce about ten 
thousand points in the parameter space. The details of 
our MCMC code can be found in our earlier paper {30| . 

The marginalized one-dimensional probability distri- 
bution function (PDF) for the the decay rate F is shown 
in Fig. [U and the contour map of Q^o vs. F is plotted in 
Fig. [2l We find that the data favors a cold dark matter 
density around Q^o — 0.186, and it is consistent with 
a stable (not decaying) dark matter. The limit on the 
decay width F are 



and 



F < 0.8 X W-^"s~'^,68.3%C.L. 



F < 4.5 X lO-^°s-\95.5%C.h. 



(13) 



(14) 



H{a')a'' 



respectively. Thus, the decay-to-neutrino rate of the 
DCDM must satisfy > 0.7 x 10^ Gyr at 95.5% C.L. 
This is a much tighter constraint than previous ones [3l| . 



3 




r [io-=°s-'] 

FIG. 1: The PDF of the dark matter decay-to-neutrino rate 
obtained by using the SN la and WMAP data. 
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tory of the Universe. Using a high-quality SN la data 
set selected from several resent observations (including 
GoldOe, SNLS and ESSENCE) at low and intermediate 
redshift, and the position of the first peak of the CMB 
angular spectrum given by the WMAP three- year data at 
high redshift, we obtain a model- independent constraint 
on the decay width F by employing the MCMC tech- 
nique. We find r < 4.5 X lO-^Og-i at 95.5% C.L. i.e., 
> 0.7 X 10^ Gyr. One could put further constraint 
on the decay life time of the particle if a is known. For 
example, if a = 0.05, then the decay life time F"^ > 14 
Gyr, i.e. greater than the current age of the Universe. 

Recently the terrestrial atmospheric neutrino spectra 
observed by several neutrino detectors such as Frejus 
[H, AMANDA [3l,IM[33 and Super-Kamiokande jsl 
has been used to constrain the annihilation [s^, IH, [S^ 
or decay [i^ properties of CDM. It is assumed that 
the annihilations or decays could produce very energetic 
neutrinos, which would produce a peak in the observed 
atmospheric neutrino background. These studies gen- 
erally yields tighter constraint on the decay width of 
the dark matter than obtained here. For example, for 
1 — lOGeV, the atmospheric neutrinos yields a con- 
straint which is several order of magnitudes stronger than 
ours. However, there are certain limitations on the ap- 
plicability of such method. If the neutrinos produced in 
the decay is of relatively low energy, below the threshold 
energy of the detector, then they would not be detected 
in atmospheric neutrino experiments. Also, if the de- 
cay process involve multiple particles and produce neu- 
trinos with a spread of energy distribution, then one may 
not always be able to identify a sharp peak in the ato- 
mospheric neutrino experiment data, although in some 
analyses, e.g. that of Ref . [40| . this effect is not obvious 
due to the large bin size used. Our method is not affected 
by these problems and is therefore a more conservative 
limit which complements their studies. 
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FIG. 2: The contour map of vs. F from SN la and 
CMB data.The la(68.3%), 2cr(95.5%) and 3cr(99.7%) confi- 
dence levels are marked by red, green and blue solid lines 
respectively. 

III. CONCLUSION 

We considered a generic decaying cold dark matter 
(DCDM) model which decays invisibly, producing rela- 
tivistic neutrinos in the process. As the produced neutri- 
nos has an equation of state of 1/3 whereas the original 
cold dark matter has an equation of state of (almost) 0, 
the cosmic expansion rate is affected. We constrain the 
decay-to-neutrino rate F by analyzing the expansion his- 
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